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The Lund dipole model DIPSY is based on BFKL evolution and satura- 
tion. It can be applied to collisions between electrons, protons, and nuclei. 
In this talk I present some recent results for exclusive final states in inelas- 
tic collisions, a method to generate final states in diffractive excitation, and 
some results for collisions with nuclei. 
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Introduction 

The Pythia MC-model is the most successful description of inelastic 
reactions in DIS and pp collisions. It does, however, need input structure 
functions determined by data, and it also uses simplified assumptions about 
correlations and diffraction (with parameters retuned at different energies). 
Our aim is not to obtain the most accurate description, but instead to un- 
derstand underlying dynamics in more detail, also at the cost of lower pre- 
cision. The results are obtained in collaboration with Christoffer Flensburg 
and Leif Lonnblad. The outline is the following: 

1. Evolution of parton densities to small x 

2. Correlations and fluctuations 

3. Diffraction 

4. Nucleus collisions 



* Talk presented at International Symposium on Multiparticle Dynamics, Kielce, 
Poland, 16-21 September 2012. 
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1. Small x evolution 

1.1. Dipole cascade models 
Mueller's dipole model 

Mueller's dipole cascade model pQ is a formulation of BFKL evolution 
in transverse coordinate space. Gluon radiation from the color charge in 
a parent quark or gluon is screened by the accompanying anticharge in 
the color dipole. This suppresses emissions at large transverse separation, 
which corresponds to the suppression of small k±_ in BFKL. For a dipole 
with transverse coordinates (x, y), the probability per unit rapidity (Y) for 
emission of a gluon at transverse position z is given by 

dV a 2 (x - y) 2 _ 3a s 

— — : = — a z- -77- -7T, with a = . 1 

dY 2vr (x-z) 2 (z-y) 2 ' ir v ' 

The dipole is split into two dipoles, which (in the large -/V c limit) emit new 
gluons independently. The result is a chain of dipoles, where the number of 
dipoles grows exponentially with Y. 

When two cascades collide, a pair of dipoles with coordinates (xj, yj) and 
(xj,yj) can interact via gluon exchange with the probability where 



J2 



fij = f(xi,yi\xj,yj) = -rf 



. / ( x i - yj) (yi - xjY 



log 



(xj -xj) 2 (yi -yj. 



(2) 



Summing over all dipoles in the cascades then reproduces the LL BFKL 
result. 

Lund cascade model DIPSY 

The Lund dipole cascade model DIPSY [21 [3] is a generalization of 
Mueller's model, which includes: 

- Important non-leading BFKL effects. The most essential ones are 
related to energy conservation and the running coupling. 

- Saturation from pomeron loops in the evolution. This is not included 
in Mueller's model or in the BK equation. 

- Confinement effects. Needed to satisfy the Froissart bound. 

- MC implementation DIPSY. This gives also fluctuations and correla- 
tions. 

- Applications to collisions between electrons, protons, and nuclei. 

At high energies several dipole pairs can interact, and in the eikonal 
approximation the amplitude is given by 

T = l-e~ F , with F = J2fij, (3) 
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where the Born amplitude F is given by summing over pairs of a dipole i 
in the projectile and j in the target. The total and elastic cross sections 
are then given by da e i/d 2 b = T 2 , da t ot/d 2 b = 2T. Multiple interactions 
give color loops, which are related to pomeron loops. In the DIPSY model 
saturation effects from color loops within the cascade evolution are also 
included, which makes the result approximately independent of the Lorentz 
frame used for the analysis. 

The DIPSY model and the MC can be applied to collisions with elec- 
trons, protons and nuclei. The coupling of a virtual photon to a qq dipole is 
determined by QED. For a proton we assume an initial wavefunction repre- 
sented by three dipoles forming an equilateral triangle. When this is evolved 
to small x-values, the result is not sensitive to the details of the initial state. 
For a nucleus the nucleons are located according to a Wood-Saxon distri- 
bution. Results for inclusive and (quasi) elastic observables were presented 

inRefs. [2 nam. 

The model can also be used for exclusive final states. BFKL describes 
results for inclusive observables, while the CCFM 0E] model, and its gen- 
eralization the LDC model [U[8], can be used for a description of exclusive 
states. Some results from Ref. [9] are shown in Fig. [1] 




VL (leading Irock] |C,eV| 



Fig. 1. Min. bias r\ distributions at 0.9 and 7 TeV (left and center). Underlying 
event at 7 TeV: N c h in transverse region vs p± for a leading particle (right). Data 
from ATLAS Q33QI]. 



2. Correlations and fluctuations 

In the MC it is also possible to calculate correlations, e.g. double parton 
distributions relevant for multiple interactions. We define the double parton 
distribution T and impact parameter profile F by the relation 

r(xi, x 2 , b; Ql Ql) = D(x u Q 2 ) D(x 2 , Ql) F{b; x 1 ,x 2 ,QlQ 2 2 ). (4) 

Some results from Ref. |12] are shown in Fig. [2] left. We note that for larger 
Q 2 hotspots develop with increased correlation at small separations b. 
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Fig. 2. Left: Two-parton impact parameter profile for proton evolution to y = at 
y/s = 15 TeV. Right: ^-distribution in pp collisions at y/s = 546 GeV and (Mx) = 
140 GeV, compared to data from UA4 [15] . 

Fluctuations are very large in QCD evolution. Taking them into account 
modifies unitarization effects. When the interaction probability fluctuates 
from event to event, the elastic amplitude is given by (T) = (l — e~ F ), which 
is not equal to 1 — e~( F \ This effect suppresses interaction for small b, and 
enhances it for larger b- values [3J. Fluctuations also give odd eccentricity 
moments, such as triangular flow in pp and AA collisions (see Refs. [13, 14J). 

3. Diffraction 

Fluctuations also cause diffractive excitations, as described in the Good- 
Walker formalism. If the projectile has an internal structure, the mass 
eigenstates ^ can differ from the eigenstates of diffraction <3? n , which have 
different eigenvalues T n . For an incoming proton the diffractive state is then 
a different mixture of the mass eigenstates, and the inclusive cross section 
for diffractive excitation is determined by the variance, (T 2 ) — (T) 2 , of the 
amplitude. Assuming that the diffractive eigenstates are represented by 
parton cascades, some results for cross sections da / dln(M^) in DIS and pp 
scattering are presented in Ref . [15] . The results do reproduce experimental 
data from HERA and the Tevatron, and we want to emphasize that the 
model is tuned only to o"t t an d o~ c \, with no new free parameter. 

We also note that for pp collisions the fluctuations are very much sup- 
pressed for central collisions, when the black limit is approached. Therefore 
diffractive excitation is largest for peripheral collisions, in a circular ring 
expanding slowly to larger radius at higher energy. 

Relation Good Walker vs triple-pomeron 

Diffractive excitation to high masses is more traditionally described in 
the triple-regge formalism. In Ref. [IB] it is, however, shown that the Good- 
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Walker and triple-pomeron formalisms are just different ways to describe the 
same phenomenon. Such a relation was also indicated in Ref. [15] . where it 
was demonstrated that the DIPSY results have the expected triple-pomeron 
form. Switching off the saturation effects in the simulation, the resulting 
bare pomeron corresponds to a single pomeron pole, with an almost con- 
stant triple-pomeron coupling. The Good-Walker formalism has then the 
advantage that the result can be calculated without extra free parameters, 
also including saturation effects, which in the triple-regge formalism are 
described by "enhanced diagrams" . 

Exclusive final states in diffractive excitation 

Diffr active excitation is fundamentally a quantum effect. For exclusive 
final states different contributions interfere destructively, and therefore there 
is no probabilistic description. Still, it is demonstrated in Ref. [17] how the 
different contributions can be calculated within the DIPSY MC, added and 
squared, to give the desired cross section. Some early applications to DIS 
andpp collisions are presented in Ref. p. 7]. As an example Fig. [2] right shows 
the pseudorapidity distribution in pp collisions at 546 GeV. We note here 
that the ^-distribution is somewhat too hard for large rj. This is due to the 
lack of quarks in the proton wavefunction, and the resulting absence of a 
forward baryon. This has to be changed in future improvements. I want to 
emphasize again that these results are based purely on fundamental QCD 
dynamics, with no free parameters beyond those tuned to the total and 
elastic cross sections. 

4. Nucleus collisions 

The DIPSY model can also be applied to nucleus collisions. Here it gives 
the full partonic picture, accounting for saturation within the cascades, cor- 
relations and fluctuations (it gives e.g. triangular flow [H]), and finite size 
effects. As example Fig. [3] left shows the shadowing effect in pPb collisions. 

For the final state in nucleus collisions the model gives a dense gluon 
soup. This state can be used as initial conditions in a subsequent hydro- 
dynamical expansion. Only adding FSR and hadronization as for hadronic 
collisions, would give too many particles, but a toy model, in which the 
gluons within 1 fm are allowed to interact and reconnect, might simulate 
the "thermalization" in the final state. As example Fig. [3] right shows the 
resulting multiplicity distribution for \rj\ < 0.5 in central collisions at RHIC 
and LHC. 

5. Summary 

The DIPSY model is based on QCD dynamics for small x evolution plus 
saturation. It is an attempt to understand basic dynamics, not to obtain 
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Fig. 3. Left: Shadowing effect R = a p A/Aa pp for pPb collisions. Right: Charged 
multiplicity from toy model "thcrmalization" , in |ry| < 0.5 for central collisions at 
RHIC (AuAu) and LHC (PbPb). 



optimal precision. 

- It works well for inclusive observables, 

- gives a fair description of exclusive final states, 

- describes diffraction with no extra free parameter (including exclusive 
diffractive states), 

- includes correlations and fluctuations (goes beyond the meanfield ap- 
proximation in the BK eq.), 

- is also applicable to nucleus collisions, where it accounts for saturation 
within the cascades, correlations, fluctuations, and finite size effects, and can 
give initial conditions for hydro expansion. 
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